Construction of a High-Density Genetic Linkage Map for the Mapping of QTL Associated with Growth-Related Traits in Sea Cucumber (Apostichopus japonicus)

Simple Summary Slow growth and germplasm degradation have restricted the sustainable commercial development of the sea cucumber industry. To analyze the genetic mechanism of growth traits of sea cucumbers, we constructed a high-density genetic linkage map based on single nucleotide polymorphism (SNP) molecular markers and performed a quantitative trait loci (QTL) mapping analysis. We annotated a critical candidate gene related to growth traits and explored mRNA expression levels. The results showed that the gene was significantly highly expressed during the larval developmental stages. These results can be used to genetically improve the growth traits of sea cucumbers. Abstract Genetic linkage maps have become an indispensable tool for genetics and genomics research. Sea cucumber (Apostichopus japonicus), which is an economically important mariculture species in Asia, is an edible echinoderm with medicinal properties. In this study, the first SNP-based high-density genetic linkage map was constructed by sequencing 132 A. japonicus individuals (2 parents and 130 offspring) according to a genotyping-by-sequencing (GBS) method. The consensus map was 3181.54 cM long, with an average genetic distance of 0.52 cM. A total of 6144 SNPs were assigned to 22 linkage groups (LGs). A Pearson analysis and QTL mapping revealed the correlations among body weight, body length, and papillae number. An important growth-related candidate gene, protein still life, isoforms C/SIF type 2 (sif), was identified in LG18. The gene was significantly highly expressed during the larval developmental stages. Its encoded protein reportedly functions as a guanine nucleotide exchange factor. These results would facilitate the genetic analysis of growth traits and provide valuable genomic resources for the selection and breeding of new varieties of sea cucumbers with excellent production traits.


Introduction
Sea cucumber (Apostichopus japonicus) is an important edible aquaculture species in East Asia because of its considerable medicinal value. In 2020, more than 171,000 tons of sea cucumber were produced from a farming area of 246,745 hectares. The cultivation tive genomic regions and target low-copy regions 2-to 3-times more efficiently, is the key to a successful GBS analysis. Compared with the RAD method, the GBS method is less complicated. The single-well digestion of genomic DNA and ligation of adapters result in decreased sample processing, fewer DNA purification steps, and the elimination of a fragment-size selection step [13]. In the present study, we used an F 1 population as the experimental material to construct the first high-density A. japonicus genetic linkage map with 6144 SNPs by applying GBS technology. Body weight, body length, and papillae number were determined for each individual and used for quantitative trait locus (QTL) mapping. Eleven growth-related candidate genes were identified. This map is a useful resource for future research on mining of information on functional genes associated with economically important traits, as well as for implementing molecular MAS breeding of A. japonicus.

Mapping Population and DNA Extraction
A mapping population was established by mating a pair of unrelated and highly polymorphic male (MP) and female (FP) sea cucumber specimens (with a wet weight of 96 g and 235 g, body length of 17.2 cm and 24.3 cm, respectively), which were collected from a wild population in Rushan, Shandong province, China. After the fertilization, the body wall and muscle tissues of the parents were sampled and preserved in ethanol. A total of 130 randomly selected 1-year-old offspring and the parents were used for genetic mapping and QTL detection. Growth traits, including wet body weight, body length, and papillae number, were recorded for all individuals. The body wall and muscle tissues were collected from all specimens and stored at −80 • C. Genomic DNA was extracted from the frozen samples according to a traditional phenol-chloroform DNA extraction protocol [14]. The quality of the genomic DNA samples was evaluated via 1% agarose gel electrophoresis and the concentrations were determined using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). High-quality genomic DNA (i.e., in terms of purity and integrity; OD 260 /OD 280 = 1.8-2.0 and OD 260 /OD 230 = 1.8-2.0) with a concentration greater than 50 ng/µL was used to construct the GBS libraries.

GBS Library Construction and Sequencing
To ensure that we selected an appropriate restriction enzyme, which is critical for a successful GBS run, we conducted a digital enzyme cut analysis first. The GBS libraries were constructed for the 2 parents and 130 offspring according to a previously described standard protocol [13]. Briefly, we digested 0.1-1 µg genomic DNA using MseI (cut-site: TTAA), HaeIII (cut-site: GGCC), and MspI (cut-site: GGCC), which are insensitive to Dam, Dcm, and CpG methylation. Common and barcoded adapters with sticky ends were ligated to the digested fragments before the PCR amplification and Illumina sequencing. The primer sequences were as follow: Primer 1: 5 -AATGATACGGCGACCACCGAGATCTACACTCTTT CCCTACACGACGCTCTTCCGATCT-3 and Primer 2: 5 -CAAGCAGAAGACGGCATACG AGATCGGTCTCGGCATTCCTGCTGAACCGCTCTT CCG ATCT-3 [15]. The Illumina HiSeq 2500 sequencing platform was used for the paired-end sequencing of the GBS libraries. [5]. The raw read data were archived at the NCBI Sequence Read Archive (SRA) under Accession Number PRJNA777474. The raw reads for each sample were filtered during a quality control step to eliminate reads containing adapters and/or ambiguous bases (N) as well as low-quality reads (>50% bases harboring Phred quality score ≤ 20) [4,16,17].

SNP Calling and Genotyping
The retained clean reads for each individual were aligned to the A. japonicus genome [5] using BWA (version 0.1.22) [17]. The reads that were aligned to only one position of the reference genome were selected for subsequent analyses. The SAMtools (version 1.0) program was used for detecting SNPs and for filtering the BAM files [18]. To decrease the number of false positives, the SNP base support number for the parents and offspring was required to be no less than 4 and 1, respectively. To minimize the number of false positives caused by repeated regions, the SNP base support number was required to be no more than 1000. Additionally, the genotype quality threshold was set at ≥10. The genotype detection results for the FP and MP were used to identify the SNPs between the parents. The sites lacking parental information were filtered out, and those that were homozygous and polymorphic between the parents were retained. These filtered SNPs were categorized as follows: 'hk × hk' (markers in both parents), 'lm × ll' (markers in the MP), or 'nn × np' (markers in the FP). The three types of markers were further screened by checking for abnormal bases and by filtering for genotyping completeness.

Linkage Map Construction
Sex-specific linkage maps were constructed using the cross-pollination algorithm of JoinMap (version 4.0) [19]. A logarithm of odds (LOD) threshold of 3.0 was applied to assign markers to linkage groups (LGs). Recombination rates were calculated according to the regression mapping algorithm and converted into map distances (i.e., centimorgans (cM)) using the Kosambi mapping function. A consensus genetic linkage map was constructed using MergeMap [20] by integrating the sex-specific maps on the basis of shared markers. All genetic linkage maps were drawn using Perl SVG.

QTL Mapping and Potential Candidate Gene Identification
The phenotypic data (body weight, body length, and papillae number) of the 130 offspring were preprocessed using R (version 4.0.2). The growth-related QTL were mapped using the 'interval mapping' and 'restricted multiple QTL model mapping' algorithms of MapQTL (version 6.0) [21]. The LOD significance thresholds for each trait were calculated according to a permutation test (1000 permutations) at a significance level of p < 0.05. The QTL confidence interval was calculated as previously described [22]. The phenotypic variance explained (PVE) for each QTL was calculated on the basis of the offspring population variance. Furthermore, the QTL results and A. japonicus genome sequence [5] were used to identify SNP loci related to growth and for screening candidate growth-related genes.

Expression and Growth Association Analyses of a Potential Candidate Gene
A quantitative real-time polymerase chain reaction (qRT-PCR) analysis was conducted to explore the expression characteristics of a potential candidate gene during the following A. japonicus larval developmental stages: fertilized egg, blastula, gastrula, early auricularia larva, mid-auricularia larva, late-auricularia larva, doliolaria larva, pentactula larva, and juvenile [23,24]. Samples were treated with RNAwait (Solarbio, Beijing, China) for the subsequent RNA extraction [25]. Three tissues (body wall, respiratory tree, and intestine) were collected from three randomly selected individuals for the gene expression analysis. Total RNA was extracted from the collected tissues using a TaKaRa MiniBEST Universal RNA Extraction Kit (TaKaRa, Dalian, China) following the manufacturer's instructions. The quality of the RNA samples was evaluated by 1.5% agarose gel electrophoresis and the concentrations were determined using the NanoDrop 2000 spectrophotometer. Highquality RNA (i.e., in terms of purity and integrity; OD 260 /OD 280 = 1.8-2.0 and OD 260 /OD 230 = 1.8-2.0) with a concentration greater than 50 ng/µL was used for the subsequent analysis. The reverse transcription reactions and qRT-PCR analyses were performed using the Prime Script™ RT reagent Kit (TaKaRa, Dalian, China) and the SYBR ® Premix ExTaq™ system TaKaRa, Dalian, China). The reference gene cytochrome b served as the internal control. Three replicates per sample were analyzed. Relative expression levels were calculated using the 2 −∆∆Ct method [26]. All the statistical analyses were performed using GraphPad Prism (version 9.0.0.121, GraphPad Software, San Diego, CA, USA).

Genotyping-by-Sequencing and SNP Discovery
The paired-end sequencing of the GBS libraries using the Illumina HiSeq 2500 platform generated 42.61 Gb raw data for 130 offspring (327,790 bp per individual) and the FP ). The quality of the sequencing data was high (Q20 and Q30 >90%) and the GC content was normal (>36%). Thus, the data were suitable for the subsequent analyses (Table S1). The use of MseI, HaeIII, and MspI resulted in 3,532,902 bp (FP), 2,740,669 bp (MP), and 176,567,947 bp (offspring) clean reads (Table S2).
Single nucleotide polymorphisms in the analyzed population were detected using the filtered BAM files. A total of 75,996 SNPs were detected in the FP, of which 2770 and 73,226 were homozygous and heterozygous SNPs, respectively. Additionally, 45,490 SNPs were detected in the MP, of which 16,403 and 29,087 were homozygous and heterozygous SNPs, respectively. A total of 48,675 polymorphic loci between the parents were detected. After filtering for genotyping completeness, the 'lm × ll' and 'hk × hk' categories comprised 2391 and 1521 SNPs, respectively (75% coverage), whereas the 'nn × np' category consisted of 1864 SNPs (66% coverage).

Genetic Linkage Map Construction
Both sex-specific linkage maps that were constructed included 22 LGs. The female map contained 3120 SNPs and spanned 2865.97 cM, with an average genetic distance of 0.92 cM. The genetic distance of each LG ranged from 20.58 cM (LG6) to 244.90 cM (LG1), and the number of markers in each LG ranged from 72 to 311, with the exception of LG6, which contained 20 SNPs (Tables 1 and S3). The male map comprised 3814 SNPs and spanned 1746.60 cM, with an average genetic distance of 0.46 cM. The genetic distance of each LG ranged from 41.00 cM (LG21) to 116.99 cM (LG1), and the number of markers in each LG ranged from 115 to 265, with the exception of LG20 and LG21, which contained 44 and 36 SNPs, respectively (Tables 1 and S3). The marker gap analysis revealed that 98.0% of the marker gaps in the male map were less than 5 cM and no marker gap was longer than 20 cM. In contrast, 95.8% of the marker gaps in the female map were less than 5 cM and one marker gap was longer than 20 cM (Table 2). These results indicated that the markers were distributed more evenly in the male map than in the female map. The sex-specific maps were integrated into a consensus map ( Figure 1) with 22 LGs, which contained 6144 SNPs and spanned 3181.54 cM, with an average genetic distance of 0.52 cM. The genetic distance of each LG ranged from 78.95 cM (LG20) to 254.19 cM (LG1), and the number of markers in each LG ranged from 99 to 522. The marker gap analysis indicated that 98.5% of the marker gaps were less than 5 cM and no marker gap was longer than 20 cM. Thus, the markers were distributed more evenly in the consensus map than in the sex-specific maps (Tables 1 and 2). The sex-specific maps were integrated into a consensus map ( Figure 1) with 22 LGs, which contained 6144 SNPs and spanned 3181.54 cM, with an average genetic distance of 0.52 cM. The genetic distance of each LG ranged from 78.95 cM (LG20) to 254.19 cM (LG1), and the number of markers in each LG ranged from 99 to 522. The marker gap analysis indicated that 98.5% of the marker gaps were less than 5 cM and no marker gap was longer than 20 cM. Thus, the markers were distributed more evenly in the consensus map than in the sex-specific maps (Tables 1 and 2).

QTL Mapping and Growth Association Analysis
A QTL mapping analysis was performed for A. japonicus growth-related traits (body weight, body length, and papillae number). The frequency distribution of phenotypic data is shown in Figure 2. The Pearson correlation analysis revealed that there were significant correlations among these three traits ( Figure 3). The QTL mapping analysis involving the consensus map indicated that the QTLs for these three traits were all located close together in the same LGs (LG5, LG18, and LG22), which reflected the correlations among body weight, body length, and papillae number (Tables 3 and S4). A total of 22 significant QTLs were detected with the LOD threshold set as 3.0. The most significant QTL related to body weight was located at 65.19 cM, presenting the highest LOD value of 5.06, explaining the highest percentage of the phenotypic variation (PVE) of 16.5 detected on LG18; that related to body length was located at LG5 (LOD 4.47, PVE 14.8%); and for papillae number was at LG18 (LOD 4.13, PVE 13.8%).

QTL Mapping and Growth Association Analysis
A QTL mapping analysis was performed for A. japonicus growth-related traits (body weight, body length, and papillae number). The frequency distribution of phenotypic data is shown in Figure 2. The Pearson correlation analysis revealed that there were significant correlations among these three traits (Figure 3). The QTL mapping analysis involving the consensus map indicated that the QTLs for these three traits were all located close together in the same LGs (LG5, LG18, and LG22), which reflected the correlations among body weight, body length, and papillae number (Tables 3 and S4). A total of 22 significant QTLs were detected with the LOD threshold set as 3.0. The most significant QTL related to body weight was located at 65.19 cM, presenting the highest LOD value of 5.06, explaining the highest percentage of the phenotypic variation (PVE) of 16.5 detected on LG18; that related to body length was located at LG5 (LOD 4.47, PVE 14.8%); and for papillae number was at LG18 (LOD 4.13, PVE 13.8%).       Eleven growth-related genes were identified from QTL mapping and the high-density genetic linkage map based on the A. japonicus genome annotation (Tables 4 and S5). Among them, 8 genes were located in the genetic linkage region, including caveolin-1, homogentisate 1,2-dioxygenase (HGD), protein still life, isoforms C/SIF type 2 (sif ), IgGFc-binding protein (FCGBP), nerve growth factor-like protein precursor (NGF), nucleoside diphosphate kinase 7, hyalin-like, neural cell adhesion molecule 1-like (NCAM), while 3 were in intergenic regions. The gene annotated as cyclic AMP-responsive element-binding protein 3-like protein 1 (CREB3L1), 150 bp upstream of the marker lm2694, was reported to be involved in growth hormone synthesis, secretion, and action. Transposon-derived Buster3 transposase-like protein (ZBED8) was 68 bp upstream of the marker lm575, and transcription factor 4 isoform X3 was located at 1239 bp upstream of marker np845. Among them, we noticed that the markers lm863 and lm864, which were both annotated to the sif gene on scaffold 157, had the highest PVE in body weight and papillae number, which means that they explain the highest percentage of the phenotypic variation, which suggested that sif may be a candidate growth-related gene. This gene encodes a guanine nucleotide exchange factor that is mainly related to the differentiation of nerve synapses.

Expression and Growth Association Analyses of a Potential Candidate Gene
To clarify the function of the sif gene, we performed qRT-PCR analysis to examine its level of expression at different larval developmental stages as well as in different tissues of adult specimens. The qRT-PCR primers are listed in Table S6. The qRT-PCR data indicated there were significant differences in expression among the examined larval developmental stages (p < 0.01). The relative expression of sif was highest in the mid-auricularia larval stage (6-times higher than the corresponding level in the fertilized egg stage) (Figure 4a). However, no significant differences were detected among the body wall, respiratory tree and intestine of adults (Figure 4b).  of adult specimens. The qRT-PCR primers are listed in Table S6. The qRT-PCR data in cated there were significant differences in expression among the examined larval dev opmental stages (p < 0.01). The relative expression of sif was highest in the mid-auricula larval stage (6-times higher than the corresponding level in the fertilized egg stage) (F ure 4a). However, no significant differences were detected among the body wall, resp tory tree and intestine of adults (Figure 4b).

Discussion
Genetic linkage maps are a crucial resource for genetics and genomics research this study, we analyzed 130 F1 offspring and their parents using GBS technology to c struct the first high-density, SNP-based genetic linkage map of sea cucumber (A. jap cus). The map was divided into 22 LGs, which is consistent with the results of an ear investigation regarding A. japonicus chromosomes [1]. The established genetic link map provides a foundation for the mapping QTL associated with important quantitat and economically valuable traits as well as for inferring the genes responsible for th traits in A. japonicus.
A linkage map with a small marker interval (<2 cM), high sequence coverage, a uniform marker distribution may be called a high-density genetic map; such a map u ally includes thousands of polymorphic markers [7,27]. The construction of a high-den genetic map requires a segregating population derived from genetically diverse pare with considerable number of polymorphisms. Commonly used populations in genom research include F1, F2, backcross, doubled haploid, and recombinant inbred lines [ Because of the long generation time for sea cucumber, we selected the F1 population this study. More specifically, a pair of unrelated sea cucumbers with large difference weight were selected as the parents for the hybridization that produced the F1 populati In the past few years, microsatellites and AFLPs have been commonly used to constr genetic maps, most of which were low-density maps lacking sufficient and/or accur information regarding the number and locations of genes or QTL. For example

Discussion
Genetic linkage maps are a crucial resource for genetics and genomics research. In this study, we analyzed 130 F 1 offspring and their parents using GBS technology to construct the first high-density, SNP-based genetic linkage map of sea cucumber (A. japonicus). The map was divided into 22 LGs, which is consistent with the results of an earlier investigation regarding A. japonicus chromosomes [1]. The established genetic linkage map provides a foundation for the mapping QTL associated with important quantitative and economically valuable traits as well as for inferring the genes responsible for these traits in A. japonicus.
A linkage map with a small marker interval (<2 cM), high sequence coverage, and uniform marker distribution may be called a high-density genetic map; such a map usually includes thousands of polymorphic markers [7,27]. The construction of a high-density genetic map requires a segregating population derived from genetically diverse parents with considerable number of polymorphisms. Commonly used populations in genomics research include F 1 , F 2 , backcross, doubled haploid, and recombinant inbred lines [28]. Because of the long generation time for sea cucumber, we selected the F 1 population for this study. More specifically, a pair of unrelated sea cucumbers with large differences in weight were selected as the parents for the hybridization that produced the F 1 population. In the past few years, microsatellites and AFLPs have been commonly used to construct genetic maps, most of which were low-density maps lacking sufficient and/or accurate information regarding the number and locations of genes or QTL. For example, a previously published A. japonicus genetic map constructed according to an F 1 pseudo-testcross mapping strategy contained fewer than 200 sex-specific markers, with a marker density exceeding 7.0 cM [3]. Another A. japonicus genetic map, which was constructed on the basis of microsatellites and SNPs, consisted of approximately 200 markers, with a marker density of 7.0 cM [2]. Compared with traditional molecular markers, SNPs can be typed on a larger scale and they are generally more abundant in genomes, enabling analyses of all genomic regions [27]. To reveal sufficient sea cucumber genetic variation, we adopted GBS technology for genotyping. Compared with the RAD method, the GBS process is simpler because the single-well digestion of genomic DNA and adapter ligations help to minimize the required sample processing, resulting in relatively few DNA purification steps and no need for selecting fragment sizes. RAD-seq contains the step of ultrasonic random interruption, thus, the cost is higher than GBS, and RAD-seq is suitable for the study without reference genome [13]. By using the GBS approach, we constructed an A. japonicus genetic map 3181.54 cM in length, with 6144 SNPs and a map resolution of 0.52 cM. The number of markers and the resolution were substantially greater for our A. japonicus genetic map than for earlier maps constructed using more traditional markers.
In this study, the genetic distance was greater in the female map than in the male map, with a distance ratio of 1.64:1 (F:M), which was consistent with the results of previous investigations [2,4]. Regarding the species examined in earlier related studies, there were sex differences in the recombination rates of autosomes. As explained by Trivers [29], the recombination rate should be lower for males than for females because they have undergone stronger sexual selection; the meiosis in females occurs only during fertilization, and the chances of female haploid selection are relatively low. In the current study, an analysis of marker gaps detected 3738 gaps that were less than 5 cM in the male map (98.0% of all gaps), with no gaps longer than 20 cM. In contrast, there were 2978 gaps that were less than 5 cM in the female map (95.8% of all gaps), and one gap was longer than 20 cM. Additionally, the female map contained 99 gaps ranging from 5 to 10 cM and 21 gaps ranging from 10 to 20 cM, whereas the male map included only 43 gaps that were 5-10 cM and 12 gaps that were 10-20 cM. These findings imply that the male map was more uniform than the female map. A consensus map was constructed on the basis of the common markers between the MP and FP. The marker gap analysis of the consensus map indicated that 98.5% of the gaps were less than 5 cM, with no gaps greater than 20 cM. Hence, the consensus map was more uniform than the sex-specific maps.
The high-density genetic linkage map constructed enabled the mapping of QTL related to A. japonicus growth character. Growth traits are affected by multiple genes and are influenced by other factors. Moreover, they are among the most important traits affecting the economic value of A. japonicus. The genetic mechanism underlying growth has long been a topic of considerable interest among breeders. However, most research on sea cucumbers focused on stress-resistance traits, body color, and regeneration. In our study, we mapped the QTL associated with growth traits. The papillae number, as an important measure of the quality of sea cucumbers, affects the market price of sea cucumbers to a large extent [30], and is also an important breeding trait related to growth. Thus, we designated body weight, body length, and papillae number as three indicators to measure growth traits. Our analyses revealed the proximity of three growth-related QTL in LG5, LG18, and LG22. The results of Pearson's association analysis and QTL mapping reflected the correlations among body length, body weight, and papillae number. Therefore, we may be able to improve the overall growth of A. japonicus by manipulating one of these traits. We annotated 11 candidate genes related to growth, which were mainly related to nerves, signal transduction, material transport, cell proliferation, and apoptosis. CREB3L1, is a membrane-bound transcription factor of the CREB/activating transcription factor (ATF) family, expressed in the glial cells, included in growth hormone synthesis/secretion [31], and is reported as may play an important role in limiting virus spread by inhibiting proliferation of virus-infected cells [32]. Caveolin-1 is an integral membrane protein soluble in multiple cellular compartments, and is reported to have functions in membrane trafficking, lipid transport, and signal transduction [33]. Nerve growth factor-like protein may act as a growth factor and participate in the paracrine regulation of prostate growth [34] and the regulation of spermatogenesis [35]. The neural cell adhesion molecule 1-like is one of the most classical cell adhesion molecules in vertebrates, and also is involved in regulating the trafficking of the neurotransmitter receptor, receptor-mediated signaling and behavior [36]. Protein still life, isoforms C/SIF type 2 gene, encodes the SIF protein, which is localized in the submembrane region. The SIF protein is believed to function as a guanine nucleotide exchange factor that controls synaptic differentiation through its effects on the organization of the actin cytoskeleton, possibly via the activation of Rho-like GTPases [37]. Moreover, it may mediate developmental processes. In our study, the sif gene was detected by the two markers (lm863 and lm864) with the highest PVE for body weight and papillae number, and it was located within the genetic linkage region. Therefore, we consider it to be an important growth-related candidate gene, and we carried out a preliminary exploration of the function of the sif gene. We found that sif was expressed during different A. japonicus larval developmental stages, especially in the mid-auricularia larval stage, when it reached its highest level of expression. Further investigations are needed to elucidate the precise regulatory effects of the SIF protein on A. japonicus embryonic development and growth.

Conclusions
In conclusion, we constructed the first SNP-based high-density genetic linkage map of sea cucumber (A. japonicus). The map (3181.54 cM) comprised 6144 SNPs distributed in 22 LGs, with an average genetic distance of 0.52 cM. The mapping of QTL for three growth-related traits (body weight, body length, and papillae number) suggested that the examined traits are correlated. Further, we identified 11 growth-related genes, among them the sif gene which we consider as an important candidate gene influencing sea cucumber growth. We preliminarily explored the spatiotemporal expression characteristics of this gene. Our research provides an important theoretical basis for the genetic analysis of growth traits and valuable genomic resources for the selection and genetic marker-assisted breeding of new varieties of sea cucumbers with excellent traits.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/biology11010050/s1, Table S1: Information about samples for completed GBS sequencing; Table S2: Results about samples for enzyme digestion; Table S3: Information about the sex-specific genetic linkage map of Apostichopus japonicus; Table S4: Information about the QTL of growth traits; Table S5: Detailed information about the annotation of candidate genes related to growth; Table S6: Primers for qRT-PCR.